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A sizeable difference in the differential production cross section of top and antitop quarks, re- 
spectively, is predicted for hadronically produced heavy quarks. It is of order Os and arises from the 
interference between charge odd and even amplitudes respectively. For the TEVATRON it amounts 
up to 15% for the differential distribution in suitable chosen kinematical regions. The resulting inte- 
grated forward-backward asymmetry of 4-5% could be measured in the next round of experiments. 
At the LHC the asymmetry can be studied by selecting appropriately chosen kinematical regions. 
' Furthermore, a slight preference at LHC for centrally produced antitop is predicted, with top quarks 

^>«^ I more abundant at large positive and negative rapidities. 
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I. INTRODUCTION 



Heavy flavor production at hadron colliders is one of the most active fields of current theoretical and experimental 
studies. Large event rates, combined with improved experimental techniques, allow for detailed investigations of the 
properties of heavy quarks and their production mechanism at the same time. While charm production with a quark 
mass around 1.5 GeV is barely accessible to perturbative QCD calculations, bottom and a forteriori top production 
should be well described by this approach. 

Theoretical and experimental results |jl|,|| for the cross section of hadronic top production are well consistent with 
00 this expectation. Obviously, in view of the large QCD coupling, the inclusion of higher order QCD corrections in these 
calculations is mandatory for a successful comparison. Recent studies have, to a large extent, concentrated on the 
■ predictions of the total cross section and a few selected one particle inclusive distributions. In this paper a different 
CIh issue of heavy flavor production is investigated, namely the charge asymmetry, which is sensitive toward a specific 
^ ^ subclass of virtual and real radiative corrections. 

^ ' Evaluated in Born approximation the lowest order processes relevant for heavy flavor production 



q + q^Q + Q, (1) 

g + g~^Q + Q, (2) 

do not discriminate between the final quark and antiquark, thus predicting identical differential distributions also 
for the hadronic production process. However, radiative corrections involving cither virtual or real gluon emission 
lead to a sizeable difference between the differential quark and antiquark production process and hence to a charge 
asymmetry which could be well accessible experimentally. 

This asymmetry has its origin in two different reactions: radiative corrections to quark- antiquark fusion (Fig. ^ and 
heavy flavor production involving interference terms of different amplitudes contributing to gluon-quark scattering 
(Fig. I) 

g + q^Q + Q + q , (3) 

a reaction intrinsically of order a^. Gluon fusion remains of course charge symmetric. In both reactions (||) and ^ 
the asymmetry can be traced to the interference between amplitudes which are relatively odd under the exchange of 
Q and Q. In fact, as shown below in detail, the asymmetry can be understood in analogy to the corresponding one in 
QED reactions and is proportional to the color factor rf^^^. In contrast, the non-Abelian contributions, in particular 
those involving the triple gluon coupling, lead to symmetric pieces in the differential cross section. Event generators 
which do not include the full next-to-leading matrix elements cannot predict the asymmetry. 
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FIG. 1. Origin of the QCD charge asymmetry in hadroproduction of heavy quarks: interference of final-state (a) with 
initial-state (b) gluon bremsstrahlung plus interference of the box (c) with the Born diagram (d). Only representative diagrams 
are shown. 
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FIG. 2. Origin of the QCD charge asymmetry in hadroproduction of heavy quarks through flavor excitation. 



Let us briefly discuss a few important aspect of this calculation. The box amplitude for qq — > QQ is ultraviolet 
finite and the asymmetric contribution to the cross section of order is therefore not affected by renormalization, an 
obvious consequence of the symmetry of the lowest order reaction. The same line of reasoning explains the absence 
of initial state coUinear singularities in the limit —^ which would have to be absorbed into the (symmetric) 
lowest order cross section. Infrared singularities require a more careful treatment. They are absent in the asymmetric 
piece of the process in eq. (||). However, real and virtual radiation (Fig. |l]), if considered separately, exhibit infrared 
divergences, which compensate in the sum, corresponding to the inclusive production cross section. 

The charge asymmetry in the partonic reactions and (^) implies for example a forward-backward asymmetry 
of heavy flavor production in proton-antiproton collisions. In particular, it leads to a sizeable forward-backward 
asymmetry for top production which is dominated by reaction (|^), and can, furthermore, be scrutinized by studying 
tt production at fixed longitudinal momenta and at various partonic energies s. However, the charge asymmetry can 
also be observed in proton-proton collisions at high energies. In this case one has to reconstruct the tt restframe and 
select kinematic regions, which are dominated by qq annihilation or flavor excitation gq ttX. Alternatively, one 
may also study the difference in the one-particle inclusive rapidity distribution of top versus antitop, which again 
integrates to zero. 

The analysis of these effects allows to improve our understanding of the QCD production mechanism. At the same 
time it is important for the analysis of single top production through Wb fusion. This reaction is charge asymmetric 



2 



as a consequence of weak interactions. Although the final states in single top production and hadronic tt production 
are different and should in principle be distinguishable, it is nevertheless mandatory to control the charge asymmetry 
from both sources. 

The presence of charge asymmetric contributions in the flavor excitation reaction has also been noticed in j^-Q] for 
b quark production, however without any quantitative statement. The charge asymmetry was also investigated in 
In this work only real gluon emission was considered. To arrive at a finite result, a gluon energy infrared cut, Ecut 
had to be introduced which leads to arbitrary large results with a pronounced dependence on Ecut and a different 
sign of the asymmetry compared to our inclusive calculation. 

The outline of this paper is as follows. The technical aspects of the calculations will be presented in section ||. 
The asymmetric pieces of real and virtual corrections will be given, together with the numerical evaluation, the 
compensation of infrared singularities from real and virtual radiation, and the results at the partonic lev el. The 
implications for hadronic collisions, proton-proton as well as proton-antiproton, will be studied in section III, A brief 
summary and our conclusions will be given in section IV. 

The main thrust of this paper is toward the study of top quarks. Nevertheless all the results are in principle 
applicable to bottom and charm quarks. In practice, of course kinematical regions have to be selected, which are 
dominated by qq annihilation or flavor excitation. For b quarks, furthermore, the dilution effects of mixing must be 
included. We will comment on these points in more detail below. 



II. AMPLITUDES AND PARTONIC CROSS SECTION 



As we shall see below, the dominant contribution to the charge asymmetry originates from qq annihilation, namely 
from the asymmetric piece in the interference between the Born amplitude for qq QQ (Fig. |^d) and the one loop 
corrections to this reaction (Fig. |]c) , which must be combined with the interference term between initial state and final 
state radiation (Fig. |l|a.p3). The corresponding contribution to the rate is conveniently expressed by the absorptive 
contributions (cuts) of the diagrams depicted in Fig. 




(e) 



FIG. 3. Cut diagrams. 

However, only Fig. ||a plus the crossed Fig. ||b are relevant for the charge asymmetric piece. Fig. ||c, ^d and ||e, on 
the other hand lead to a symmetric contribution only. This can be seen as follows: the color factors corresponding to 
Fig. 0a and 0b (after averaging over initial and summing over final states) respectively are given by 
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^» - ^^-^ (tTt) (tTt) - T^'-^-'' " • 

where Nc = 3, /^^^ = 24 and d^^^ = 40/3. Without color factors the contributions to the cross section from Fig. ^ 
and l^b are related by 

dcr3a{Q,Q) = -ddsbiQ.Q) , (5) 

which holds true both for two and three particle cuts. The asymmetric piece thus originates from the d^^^ term, and 
its form is thus equivalent to the corresponding QED reaction with the replacement of the quark charges and QED 
coupling by the color factor 

The production cross section, on the other hand is obtained from the corresponding QED process through the 
replacement 

o^QEDQqQq -^NcTfCfoI , (7) 
with Tp — 1/2, Cp — 4/3. The QCD asymmetry is thus obtained from the QED results by the replacement 

Let us note in passing that the cuts through diagrams involving the triple gluon coupling (Fig. ||d, ^) lead 
to charge symmetric terms. For Fig. ^c and ^ this can be seen as follows: its contribution involves the factor 
Tr{'^°'{Q + M)j^{0 — M)} which is evidently symmetric under the exchange of Q and Q. The remainder of the 
diagram depends on Q + Q only, which leads to charge symmetry of the whole amplitude. The same line of reasoning 
applies to Fig. as far as the exchange of the initial quarks are concerned. Charge conjugation invariance then 
implies, that also Fig. ^ is symmetric under the exchange between Q and Q. These terms have to be combined with 
the charge symmetric CfCa terms from Fig. ||a and ^3 to yield a gauge invariant combination. 

Although the relevant ingredients for the charge asymmetric piece are already listed in the original QED publications 
(and the later works on hadronic heavy flavor production in the compact formulae shall be listed in the 

appendix for completeness and convenience of the rea der . In a first step virtual and soft radiation are combined, with 



a cut on the gluon energy, -E^uj (see appendix, eq. (A5)). The logarithmic divergence of this result for small -E^i^j 
is cancelled by the corresponding divergence for real radiation. A particularly compact formula for the asymmetric 
piece of the hard radiation is given in eq. ( [A3| ) . To obtain finally the asymmetric piece of the inclusive cross section 
for 

g + q-">Q + X, (9) 

the integral over the real gluon spectrum is performed numerically. 

The differential charge asymmetry at the partonic level can then be defined through 

i(cos^)^^-^--^:)"^^^-^), (10) 

7Vt(cos6i) + iVf(cos6') 

where denotes the top quark production angle in the qq restframe and N{cos9) — da/d^l{cos6). Since Nf(cos9) = 
Nt{— cos 6) as a consequence of charge conjugation symmetry, A{cos9) can also be interpreted as a forward-backward 
asymmetry of top quarks. In this case the denominator is of course given by the Born cross section for the reaction 
qq ^ QQ, (see eq. ([aT^)). In Fig. § i(cos^) is displayed for ^/§ = 400 GeV, 600 GeV and 1 TeV for mt = 175 GeV. 



For completeness we also display the result for bb production at \/I = 400 GeV with mf, = 4.6 GeV. The strong 
coupling constant is evaluated at the scale fi = Vs/2 from as{mz) = 0.118. 
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FIG. 5. Integrated charge asymmetry as a function of the partonic center of mass energy for top and bottom quark pair 
production. 

The integrated charge asymmetry 

|_ iVt(cos0>O)-iVt(cos0>O) 
iVt(cos ^ > 0) + Nt{cos 9>0) ' 

is shown in Fig. |^ as a function of With a typical value around 6 — 8.5% it should be well accessible in the next 
run of the TEVATRON. 

As mentioned already in the introduction, the asymmetric piece does not exhibit a divergence, even in the limit of 
vanishing initial quark mass; in other words, no coUinear singularities arise. The virtual plus soft radiation on one 
hand and the real hard radiation on the other contribute with opposite signs, with the former always larger than the 
later which explains the difference in sign between our result and @] . 
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FIG. 6. Asymmetric parts of the differential top quark pair production cross section from qq and qg initiated processes for 
fixed partonic center of mass energies \/I = 400 GeV and 1 TeV. 



Before moving to the application of these results by folding with the parton distribution functions let us first discuss 
the charge asymmetry in the quark-gluon induced reaction in eq. (^. The cross section for this reaction is obtained 
from the amplitudes depicted in Fig. |^. In fact its antisymmetric piece can be obtained by crossing directly from 
the reaction qq — > QQg and is given by eq. (A13). Again only the QED like piece contributes to the asymmetry, in 
contrast to those amplitudes induced by the triple gluon coupling. The inclusive cross section for quark production in 
quark-gluon collisions exhibits a coUinear divergence, the charge asymmetric piece is finite. The difference between Q 
and Q production (for fixed initial q and gluon directions) should not be confused with an asymmetry in the angular 
distribution of Q (or Q) , which is a trivial consequence of the asymmetric initial state configuration. 

The charge asymmetric pieces as a function of the scattering angle 



1 ( da{qq 
2 



QX) da{qq -> QX) 



dcos9 



dcos6 



_daf^ 
dcos 9 



(12) 



1 fda{qg ^ QX) da{qg ^ QX)\ _ daf 

2 \ dcosO dcosd J dcosO 



(13) 



are shown in Fig. ^ for a variety of partonic energies VI with nit = 175 GeV. 
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FIG. 7. Integrated charge asymmetric parts of the top quark pair production cross section from qq and qg initiated processes 
as a function of the partonic center of mass energy. 

The asymmetric contributions integrated in the forward-backward direction 

'^A^ —^dcose- —^dcose, i = qq,qg, (14) 
Jo a cos 6' J-i acosO 

are shown in Fig. |^. Clearly, at the TEVATRON energies with \/l < 2 TeV, the dominant asymmetric contribution 
comes from the qq initiated processes, even more so, since the quark-gluon luminosity in the relevant kinematic region 
is far below the one for quark- antiquark reactions. Furthermore, the difference between Q and Q production in quark- 
gluon collision does not exhibit a marked forward-backward asymmetry, which suppresses the qg induced asymmetry 
even further. 
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FIG. 8. Representative diagrams contributing to the QCD neutral current interference term. 

In addition to the pure QCD amplitudes also a mixed QCD-electroweak interference term will lead to an asymmetric 
contribution to the qq process. The QCD box diagram, Fig. ||c, can also give rise to tt in a color singlet configuration, 
which in turn interferes with tt production through the photon or Z (Fig. ^). A similar consideration applies to 
interference between initial and final state radiation. The resulting asymmetry is obtained from the QCD asymmetry 
through the following replacement 
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OiQED QtQq 



(1 



sl,){2Iq-AQqsl 



16s 



1 - 



(15) 



Another QED-electroweak term originates from the interference between the gluon-7 box andgluon-Z box respectively 
with the QCD Born amplitude (Fig. ||b). The result for this piece |^ is also given by eq. (p^). In total this leads to 
an increase of the asymmetry as given by pure QCD by a factor 1.09. This change is thus smaller than uncalculated 
higher order corrections. 

III. HADRONIC COLLISIONS 

The asymmetry can in principle be studied experimentally in the partonic restframe, as a function of s, by measuring 
the invariant mass of the tt system plus an eventually radiated gluon. It is, however, also instructive to study the 
asymmetry in the laboratory frame by folding the angular distribution with the structure functions 12 1. For proton- 
antiproton collisions it is convenient to consider the forward-backward asymmetry as function of the production angle 
in the center of mass system. The differential asymmetry for ^ = 2 TeV is shown in Fig. ^ which displays separately 
the contribution from qq and qg (plus qg) initiated reactions. The denominator includes both qq and gg initiated 
processes in lowest order. The numerator is evidently dominated by quark- antiquark annihilation as anticipated in |l^] 
as can be seen from Fig. ^. Inclusion of mixed QCD/electroweak interference term enhances the prediction by a factor 
1.09. 

At this point we have to emphasize that both numerator and denominator are evaluated in leading order (LO). The 
next-to-leading (NLO) corrections to the ti production cross section are known to be large |l^, around 30% or even 
more. In the absence of NLO corrections for the numerator we nevertheless stay with the LO approximation in both 
numerator and denominator, expecting the dominant corrections from coUinear emission to cancel. However, from a 
more conservative point of view an uncertainty of around 30% has to be assigned to the prediction for the asymmetry. 

For the total charge asymmetry at -^s = 1.8 TeV we predict 



A 



NticosO > 0) - jVt(cos6' > 0) 
A^t(cos6' > 0) + TVf (cos 61 > 0) 



4.8 



(16) 



where different choices of the structure function and different choices of the factorization and renormalization scale, 
jjL = mtl2 and ji = 2mt, have been considered and the factor 1.09 is included. An increase in the center of mass 
energy to 2 TeV leads to a slight decrease of our prediction to 4.6 — 5.5%. 



^This small term had been neglected in [|l3| 



8 



15 



=D 10 
(0 



O 
u 



5 

-5 
-10 
-15 



- 




- 


qq y 


qg (xlOO) 




-/ 

L ^-"^ 






— ►^^-i — 


" 1 1 1 1 1 1 1 1 


p p 

1 1 1 1 1 1 1 1 1 



-0.5 



0.5 1 
COS 9 



FIG. 9. Differential charge asymmetry in the proton-antiproton restframe, -^s = 2 TeV, using the CTEQ-1 structure function, 
fj, = rrit. The contributions from qq and qg (plus qg) initiated processes are shown separately. 

For illustrative purpose in Fig. |l^ the qq and qg{qg) induced contributions in the partonic restframe are also 
displayed separately in the x = xi — X2 and s plane for proton-antiproton collisions with n = \/S/2. Furthermore, 
as characteristical example the relative amount of giuon fusion as function of x — xi — X2 and s is shown in the two 
dimensional distribution of Fig. In regions of larger qq and correspondingly smaller gg induced reactions a larger 
asymmetry is expected. 

Bottom quark production at the LHC or TEVATRON is of course dominated by gluon fusion. The forward- 
backward asymmetry from qq and qg{qg) reactions is thus negligible, at least as far as the total cross section is 
concerned. However, the hh{g) final state with s sufficiently large, say above 300 GeV is again dominated by qq 
annihilation and a sizeable asymmetry is predicted in this kinematical region. Selecting for example VI > 300 GeV 
and I cos 6* I < 0.9 one predicts (^s = 2 TeV) 



A = 4.3-5.1% , 



(17) 



which should be accessible by experiment. A factor 0.96, from the QCD/electroweak interference, has been also 
included. 

Top-antitop production in proton-proton collisions at the LHC is, as a consequence of charge conjugation symmetry, 
forward-backward symmetric if the laboratory frame is chosen as the reference system. However, by selecting the 
invariant mass of the tt(+g) system and its longitudinal momentum appropriately, one can easily constrain the parton 
momenta such that a preferred direction is generated for quark- antiquark reactions. 

For some of the more extreme kinematic regions, namely large x and/or large s, a sizeable difference between top 
and antitop production can be observed at the LHC. In Fig. ^ the contributions from qq and qg [qg induced reactions 
are displayed separately. The production cross section per se, which is decisive for the possibility of measuring the 
asymmetry in these regions is displayed in Fig. In practice, only the region with s below 2 TeV will be observable, 
in particular at large x. 

The asymmetry, as displayed in Fig. ^ as a function of x — xi — X2 and s, is defined in the ti{g) restframe. 
From this it may seem that the reconstruction of both t, t and even the gluon is required for the study of the charge 
asymmetry in pp collisions. However, also the difference between the single particle inclusive distribution of t and i 
respectively may provide evidence for the charge asymmetry. This can be easily understood from Fig. 13, Production 



of tt{g) with negative x is dominated by initial q with small xi and q with large X2- The charge asymmetry implies 
that Q{Q) is preferentially emitted into the direction of q{q), Fig. p^a, The same line of reasoning is applicable 
for positive x, with Q{Q) again preferentially emitted in the direction of q{q), and the role of xi and X2 reversed. In 
total this leads to a slight preference for centrally produced antiquarks and quarks slightly dominant in the forward 
and backward direction, i.e., at large positive and negative rapidities. 
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FIG. 10. Contributions from qq and qg{qg) induced reactions to the charge asymmetry in proton-antiproton collisions, 
y/s = 2 TeV, as a function of xi — X2 ~ 2P2,{ttg) / y/s and s. Partonic restframe (CTEQ-1). 
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FIG. 11. Relative amount of gluon-gluon initiated processes as a function of xi — X2 — 2Pz{tig) / ^/s and s in lowest order, 
for i/s = 2 TeV in proton-antiproton collisions and = 14 TeV in proton-proton collisions. Contour lines go from 2% (white) 
to 10% (dark) in steeps of 2% in the former case and from from 10% (white) to 90% (dark) in steeps of 10% in the later 
(CTEQ-1). 
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FIG. 12. Contributions from qq and qg{qg) induced reactions to the charge asymmetry in proton-proton collisions, 
= 14 TeV, as a function of xi — X2 = 2P3{tig) / ^/s and s. Partonic restframe (CTEQ-1). 
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FIG. 13. Typical configuration of momenta for top and antitop production through quark annihilation in the region of large 
parallel momenta. 




FIG. 14. Rapidity distribution of charge asymmetry (a) and total cross section at Born order (b) of top quark production 
in proton-proton collisions, \/s = 14 TeV and fj, = mt. Contributions from qq fusion and flavor excitation, qg{qg), are shown 
separately. Laboratory frame (CTEQ-1). 
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The differential charge asymmetry 



dN{Q) dN(Q) 



dy 


dy 


dN{Q) 


dN{Q) 


dy 


+ dy 



(18) 



is shown in Fig. [l4| a for top quark production at the LHC {y/s = 14 TeV). As expected, a sizeable charge asymmetry 
is predicted in the region of large rapidity. It remains to be seen, if the low event rates in these extreme regions will 
permit the observation of this effect. The quark-gluon process is again negligible. 



IV. CONCLUSIONS 



The difference between the distributions of quarks versus antiquarks produced in hadronic collisions has been 
investigated. This asymmetry is particularly relevant for the production of top quarks in suitably chosen kinematical 
regions, but can also been observed for bottom quarks at large s. At the partonic level it amount up to 10% or even 
15%. In proton- antiproton collisions at the TEVATRON the integrated forward-backward asymmetry amounts close 
to 5%. At the LHC a slight preference for centrally produced antitop is predicted, with top quarks more abundant at 
large positive and negative rapidities. 



APPENDIX A: BASIC FORMULAE 



For completeness we summarize here the charge antisymmetric contributions to the heavy quark production cross 
section. The charge asymmetric piece of the hard gluon radiation process 



defined as 



is given by 



daf = - {d<j{qq -> QX) - d<j{qq ^ QX)) 



abc 



dy35 dyi5 dQ 4t:s IGiV^, yi2 (j/34 + 2m2) ^35 

X 1 — {via + vli + 2/23 + 2/24 + 2m2(y34 + 2m^ + ^12))+ y^^ 

- (1 ^ 2) - (3 ^ 4) + (1 ^ 2, 3 <-> 4) , 
with Nc = 3 and d^j,^ = 40/3. All the quantities are normalized to the partonic center of mass energy s, 

y^j = 2(pi • pj)/s , rm? = mg/s . 



(Al) 



(A2) 



(A3) 



(A4) 



The asymmetry is explicitly driven by the an tisy mmetric exchange of momenta {pi pj). 

On the other hand, soft radiation from eq. (Al) integrated in phase space up to a cut in the soft gluon energy, 
plus the virtual corrections to the Born process qq QQ contribute to the asymmetry as 



dcoaO 



(1 



„2 I A ^2 



4to ) 



4 log 



log(2w) + D{c) - D{-c) 



(A5) 



with 



(3 = Vl - 4m2 , c = /3 cos ^ , w = 



(A6) 
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and the functions B{c), coming from the box contribution, and D{c), from soft radiation, defined as 



B(c) 



1-c2-8to2 fl~C\ , „ 2x 



2ii2 1 



1 -c 



log' 



1 - c 



+ 



4c 2 - - 7m 



/32 (1 - 2m2) 



2^2 _ c2 \og{m ) + - log (to ) 



2^(l-8TO2 + 8m4) 



log' 



1"/? 
1 + /3 



4Li, - 



1 + 



TT 

T 



where 



log' 



1-2/ 



Li2 



]-Ll2[ + i«2 - 

1-2// \ y J \V 

1 



1-y 



1 - c 



i?e Li2{x') + - log2(x2) - log(a;2) log(l - x') 



^2(1 -c-2to2) 
In the hmit to ^ these functions simphfy considerably 



V=-[l-P+ V2(1-c-2to2; 



B{c) 



1 + c — c log 



1 - c 



Dfc) = loe 



1 - c 



- 2Li2 



log 



1 - c 



1 - c 



and become also free of final state collinear divergences since only integrable divergences appear at c 
The charge asymmetric contribution of the flavor excitation process 



defined as 



is given by 



da^^ = - {da{qg QX) - da(qg ^ QX)) , 



dyss dy45 dfl Ans 167V^ yis (2/34 + 2to2) 7/23 



c 

2/13 2/35 
2/12 2/25 

2 



(2/?3 + 2/14 + 2/35 + 2/45 + 2TO'(y34 + 2to2 - 2/15)) 



+4to" (^45 + 2/14) ^-(3^4) 



It is infrared finite and can be obtained just by crossing of momenta from eq. ( |A3[ ). 

For the convenience of the reader we also list the Born cross section for qq and gg fusion 



dcos9 



— ^(1 + + Am'^) 



Nc 2s 



da99-*QQ 
dcos9 

where Tp = 1/2 and Cp = 4/3. 



7r/3 / 1 Tf 
2J \Nc{l-c^) ~ 2CV 



1 + + Sm^ 



32to'' 
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